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BUCKLING  OF  ECCENTRICALLY  STIFFENED  ORTHOTROPIC 
CYLINDERS  UNDER  PURE  BENDING 

By  David  L.  Block 
Langley  Research  Center 

SUMMARY 

^he  stability  of  stiffened  orthotropic  cylinders  which  include  the  effects  of  eccentric 
stiffeners  and  are  loaded  with  pure  bending  or  any  combination  of  bending  and  compres¬ 
sion  is  investigated  analytically  by  means  of  a  small  deflection  theory.  Solutions  to  the 
buckling  equations  are  obtained  for  simple  support  boundary  conditions  by  use  of  the 
Galerkin  method.  The  pure  bending  buckling  loads  are  compared  with  existing  pure  com¬ 
pression  buckling  loads  for  three  contemporary  cylinder  configurations.  Also  shown  are 
predicted  buckling  modes  and  typical  interaction  curves  between  bending  and  compression^ 
The  results  show  that  eccentricity  effects  are  substantial  and  that  the  maximum  pure 
bending  load  may  be  as  much  as  40  percent  greater  than  the  pure  compression  load. 

INTRODUCTION 

The  effect  of  eccentricity  or  one-sidedness  of  stiffening  elements  in  determining 
buckling  loads  of  stiffened  cylinders  has  been  shown  to  be  of  great  importance.  (See 
refs.  1  to  4.)  In  reference  3,  a  theory  for  eccentrically  stiffened  orthotropic  cylinders 
employing  Donnell-type  assumptions  and  neglecting  prebuckling  deformations  was  derived 
from  energy  principles,  and  buckling  solutions  were  obtained  for  a  simply  supported  cyl¬ 
inder  under  any  combination  of  uniform  axial  load  and  lateral  pressure. 

The  purpose  of  this  paper  is  to  extend  the  analysis  of  reference  3  to  include  the 
buckling  of  an  eccentrically  stiffened  orthotropic  cylinder  subject  to  a  loading  of  pure 
bending  or  combination  bending  and  compression  (nonuntform  axial  loading).  The  equa¬ 
tions  developed  are  also  valid  for  an  unstiffened  cylindrical  orthotropic  shell.  Solutions 
to  the  buckling  equations  corresponding  to  classical  simple  support  boundary  conditions 
are  obtained  by  use  of  the  Galerkin  method  in  a  manner  similar  to  that  employed  in  refer¬ 
ence  5  in  which  an  unstiffened  isotropic  cylinder  under  pure  bending  was  investigated. 
Maximum  compressive  stress  resultant  at  buckling  due  to  bending,  hereafter  referred  to 
as  maximum  bending  loads,  are  compared  with  uniform  compression  loads  for  three 
types  of  cylinders  with  contemporary  geometry,  that  is,  ring- stiffened  corrugated  cylin¬ 
ders,  ring-  and  stringer-stiffened  isotropic  cylinders,  and  longitudinally  stiffened 


isotropic  cylinders.  Interaction  curves  between  bending  and  compression  and  the  buck¬ 
ling  mode  shapes  are  also  presented. 

SYMBOLS 

The  units  used  for  the  physical  quantities  in  this  report  are  given  both  in  the 
U.S.  Customary  Units  and  in  the  International  System  of  Units  (SI).  The  relationship 
between  these  two  systems  of  units  can  be  found  in  reference  6. 

A  cross-sectional  area  of  stiffener 

a  length  of  stiffened  cylinder  (see  fig.  1) 

an,bn,Cn  constant  coefficients  of  deflection  functions 

Djj,Dy  bending  stiffnesses  of  orthotropic  plate  in  longitudinal  and  circumferential 
directions,  respectively 

Djjy  twistii^  stiffness  of  orthotropic  plate 

d  stringer  spacing  (see  fig.  1) 

E  Yoimg's  modulus 

Ejj-jEy  extensional  stiffnesses  of  orthotropic  plate  in  longitudinal  and  circumferential 
directions,  respectively 

G  shear  modulus 

Gjjy  inplane  shear  stiffness  of  orthotropic  plate 

Iq  moment  of  inertia  of  stiffener  about  middle  surface  of  shell 

i,j,n,q,p  integers 

J  torsional  constant  for  stiffener 

I  ring  spacing  (see  fig.  1) 

m  number  of  half  waves  in  cylinder  buckle  pattern  in  longitudinal  direction 

Njj  total  compressive  stress  resultant 

Ny  maximum  compressive  stress  resultant  due  to  bending 

•^b 

critical  maximum  compressive  stress  resultant  due  to  bending  when  =  0 
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x,y,z 


Mx 


uniform  compressive  stress  resultant 
critical  compressive  stress  resultant  when  Nx^  =  0 
radius  of  cylinder  to  middle  surface  of  shell  (see  fig.  1) 
thickness  of  isotropic  cylindrical  shell  wall 

A 

effective  wall  thickness  of  stiffened  isotropic  cylinder,  * 

displacements  in  x-,  y-,  and  z-directions,  respectively,  of  a  point  in  middle 
surface  of  shell 

orthogonal  curvilinear  coordinates  with  origin  lying  in  middle  surface  of  shell 
(see  fig.  1) 

distance  from  centroid  of  stiffener  to  middle  surface  of  shell  (see  fig.  1), 
positive  if  stiffener  lies  on  external  surface  of  shell  and  negative  if  stiff¬ 
ener  lies  on  internal  surface 

Kronecker  delta,  6^  =  0  when  i  j  and  6jj  =  1  when  i  =  j 

Poisson’s  ratios  for  bending  of  orthotropic  plate  in  longitudinal  and  circum¬ 
ferential  directions,  respectively 

Poisson's  ratios  for  extension  of  orthotropic  plate  in  longitudinal  and  circum¬ 
ferential  directions,  respectively 


Subscripts: 

s,r  stringers  (longitudinal  stiffening,  parallel  to  x-axis)  and  rings  (transverse 

stiffening,  parallel  to  y-axis),  respectively 

x,y  longitudinal  and  circumferential  directions,  respectively 

A  subscript  preceded  by  a  comma  denotes  partial  differentiation  with  respect  to  the 
subscript. 

THEORETICAL  SOLUTION 

The  equations  of  equilibrium  governing  the  buckling  of  a  stiffened  cylinder  subject 
to  inplane  axial  stresses  only  are  (ref.  3): 


1  -  u  ’  d  ^(^’xx  ■  ^s^’xxx)  +  1  _  ^  t  ,\^’xy  R  j  ^xy{'^’yy  +  '^’xy)  ° 
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When  equations  (1)  for  a  stiffened  cylinder  (fig.  1)  were  developed,  the  following 
basic  assumptions  were  made.  The  stiffened  cylinder  is  composed  of  an  orthotropic  shell 
with  infinite  transverse  shear  stiffnesses  and  with  elastic  constants  as  given  in  refer¬ 
ence  7.  The  stiffening  elements  are  uniform  equally  spaced  rings  and  stringers  which 
are  assumed  to  have  elastic  properties  and  to  be  closely  spaced  so  that  their  elastic  prop¬ 
erties  can  be  averaged  over  the  stiffener  spacing.  The  usual  Donnell-type  assumptions 
are  used  to  specify  buckling  displacements  in  the  shell.  The  stiffeners  are  treated  as 
beam  elements  with  stiffener  twisting  accoimted  for  in  an  approximate  manner.  In  cases 
where  both  rings  and  stringers  lie  on  the  same  surface  of  the  shell,  the  effect  of  joints  in 
the  stiffening  framework  is  ignored. 

The  cylinder  loading  is  taken  to  be  a  combination  of  bending  and  compression,  in 
which  the  inplane  loading  prior  to  buckling  is  given  by 

Nv  =  Nx  +Nx  cos  I-  (2) 

where  Nx  is  the  uniform  compressive  stress  resultant  and  Nx^  is  the  maximum 
bending  stress  resultant;  both  are  positive  in  compression. 

The  boundary  conditions  to  be  satisfied  are  those  of  simple  support  which  are  given 
in  reference  3  as 

w  =:  0 
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Expressions  for  the  displacements  which  satisfy  these  boundary  conditions  are 


u^cosHS^  c^cosf 


oo 

v  =  sin«^  y  b„sin3J^ 

a  A  n  T? 


w  =  sin 


where  m  denotes  the  number  of  axial  half  waves  in  the  buckle  pattern. 

By  substituting  equations  (3)  into  equations  (la)  and  (lb),  the  coefficients  bj^  and 
Cj^  may  be  determined  as  functions  of  the  coefficients  and  from  equations  (Ic),  (2), 
and  (3)  the  following  equation  is  obtained: 


oo 

Q(W)  =  ^  aj^  Agg  +  A 


ri2^13  ~  •^11-^23)  ^  .  ri2^23  ~  -^13^22 

23I  ^  -  A  ^  j  A  A  -  A  ^ 

\^11^22  ^12  /  \^ll-^22  ^12  y 


sixS^cosf.O 


in  which 


Es^s 


1-MxV  ^ 


5 


where  the  reciprocal  relations  IJ-Je  =  iiJE  and  ]U.„D„  =  P„D„  (ref,  7)  have  been 

■^y  yx 

employed  to  simplify  the  equation. 

By  use  of  the  Galerkin  method,  a  set  of  equations  for  the  unknown  coefficients  a^ 
is  then 


from  which  a  homogeneous  set  of  eqxiations  is  determined  as 


where 


V  ^11^22  ■  ^12^  / 


+  Ai3 


^12^23  ~  -^13-^22^ 
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n  =  0,  1,  2,  3,  .  . 


G  =  n) 

(j  ^n) 

Equation  (6)  can  be  used  to  determine  buckling  loads  under  pure  bending  and  pure 
compression  (where  the  equation  reduces  to  the  result  given  in  ref.  3)  and  to  investigate 
interactions  between  bending  and  compression.  For  nontrivial  solutions  of  equation  (6), 
the  determinant  of  the  coefficients  must  equal  zero.  When  computing  the  buckling  load 
for  a  specified  loading,  the  determinant  is  minimized  numerically  for  integral  values 
of  m,  the  number  of  axial  half  waves  in  the  buckle  pattern.  The  size  of  the  determinant 
is  determined  by  requiring  that  the  buckling  load  converges  to  a  desired  accuracy.  For 
the  calculations  presented  herein,  the  size  of  the  determinant  was  increased  vmtil  the 
buckling  load  between  a  specified  order  determinant  and  a  determinant  of  one  less  order 
did  not  vary  by  more  than  0.5  percent.  All  the  numerical  calculations  presented  herein 
were  performed  on  an  electronic  computer. 

RESULTS  AND  DISCUSSION 

Because  of  the  large  number  of  geometric  parameters  involved  in  equation  (6),  it  is 
impractical  to  present  results  of  a  general  nature.  However,  it  is  of  value  to  present 
some  computed  results  for  cylinders  of  contemporary  proportions  in  order  to  study 
loadings  of  pure  bending  compared  with  pure  compression,  eccentricity  effects  under 
these  loadings,  and  interaction  between  bending  and  compression.  Therefore,  computa¬ 
tions  of  bending  loads  were  made  for  three  types  of  stiffened  cylinders  appropriate  for 
large -diameter  booster  interstage  structures:  ring- stiffened  corrugated  cylinders,  ring- 
and  stringer-stiffened  isotropic  cylinders,  and  longitudinally  (stringer)  stiffened  isotropic 
cylinders.  For  the  calculations,  the  material  in  the  cylindrical  shell  and  the  material  in 
the  stiffeners  were  taken  to  be  identical,  and  a  value  of  0.32  was  assigned  to  Poisson's 
ratio  for  the  stiffened  isotropic  shells.  The  dimensions  of  the  cylinders  are  given  in  fig¬ 
ure  2.  Eccentricity  effects  for  each  cylinder  configuration  were  studied  by  moving  the 
stiffeners  from  the  internal  to  the  external  surface  of  the  shell. 

The  results  of  the  computations  are  presented  in  table  I  and  in  figures  3  to  9. 
Details  and  discussion  of  the  computed  results  for  each  cylinder  configuration  follow. 

The  maximum  pure  bending  loads  for  general  instability  failure  of  ring-stiffened 
corrugated  cylinders  and  for  general  and  panel  instability  failure  of  ring-  and  stringer- 
stiffened  cylinders  are  presented  in  table  I  and  are  shown  in  figures  3,  4,  and  5,  respec¬ 
tively.  In  the  figures,  a  nondimensional  buckling  load  parameter  or  N^^/^x 
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for  corrugated  cylinders  and  or  Nx^yEt  for  ring- and  stringer-stiffened 

cylinders,  is  plotted  against  a  nondimensional  ring  spacing  l/R.  For  the  panel  insta¬ 
bility  calculations,  the  ring  spacing  I  was  taken  to  be  the  length  of  the  cylinder,  and 
thus  the  cylinder  is  a  longitudinally  (stringer)  stiffened  isotropic  cylinder.  Results  are 
presented  for  cylinders  loaded  in  pure  bending  (computed  by  using  eq.  (6))  or  loaded  in 
pure  compression  (computed  by  using  eq.  (15)  or  (17)  of  ref.  3). 

For  the  ring-stiffened  corrugated  cylinders  studied,  the  computed  results  (table  1(a) 
and  fig.  3)  indicate  that  eccentricity  effects  are  as  important  for  pure  bending  loads  as 
those  observed  for  pure  compression  (ref.  3).  For  internally  stiffened  cylinders,  the 
critical  maximum  bending  load  was  found  to  be  as  much  as  40  percent  greater  than  the 
critical  pure  compression  load.  (See  last  column  of  table  1(a).)  This  percentage  increase 
is  somewhat  greater  than  similar  results  of  25  percent  presented  in  reference  4;  how¬ 
ever,  the  corrugated  cylinders  in  reference  4  do  not  have  the  same  properties  as  the  ones 
examined  herein  and,  therefore,  the  two  results  would  not  be  expected  to  agree. 

For  the  ring-  and  stringer -stiffened  and  the  longitudinally  stiffened  isotropic  cylin¬ 
ders,  the  computed  results  (tables  1(b)  and  1(c)  and  figs.  4  and  5)  show  that  the  differences 
between  external  and  internal  stiffeners  are  not  as  large  as  those  observed  for  the  corru¬ 
gated  cylinders,  but  nevertheless,  both  effects  may  be  substantial  and  must  be  considered. 

In  figures  6  and  7,  typical  interaction  curves  between  bending  and  compression  are 
shown  for  the  corrugated  cylinder  and  for  the  ring-  and  stringer -stiffened  cylinder, 
respectively.  The  interaction  curve  for  the  internally  stiffened  corrugated  cylinder  pos¬ 
sesses  a  vertical  tangent  at  the  abscissa  as  observed  in  reference  4;  however,  for  the 
other  interaction  curves  this  behavior  was  not  as  pronounced.  Hence,  except  for  the 
internally  stiffened  corrugated  cylinder,  the  interaction  curves  presented  may  be  approx¬ 
imated  by  straight  lines. 

In  figures  8  and  9,  the  buckling  mode  shapes  for  pure  bending  and  pure  compression 
are  shown  (corresponding  to  the  end  points  of  the  interaction  curves  of  figs.  6  and  7).  The 
buckle  shapes  for  pure  bending  show  that  the  buckles  are  localized  on  the  compression 
area  of  the  cylinder,  a  result  observed  in  references  4  and  5. 

In  addition  to  these  cylinder  configurations,  calculations  were  made  for  unstiffened 
isotropic  and  orthotropic  cylinders  by  setting  the  ring  and  stringer  properties  equal  to 
zero  and  using  the  appropriate  cylindrical  shell  constants  in  equation  (6).  For  isotropic 
cylinders,  the  calculations  showed  that  the  effects  of  bending  are  negligible  when  compared 
with  compression  which  agreed  with  references  5  and  8.  The  three  cylinder  configurations 
considered  herein  were  analyzed  as  orthotropic  shells  by  using  orthotropic  constants 
defined  in  reference  3.  The  calculations  for  these  cylinders  showed  that  even  though 
eccentricity  effects  are  neglected,  there  can  be  a  considerable  increase  in  the  critical 
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compression  load  due  to  bending;  for  example,  the  corrugated  cylinder  gave  an  increase 
of  18  percent. 

Additional  information  can  be  obtained  by  considering  the  circumferential  buckling 
modes  for  pure  compression  which  are  given  in  reference  3.  The  results  presented  in 
table  I  indicate  generally  that  the  differences  between  pure  bending  and  pure  compression 
loads  may  be  considered  small  when  the  cylinder  has  an  axisymmetric  pure  compressive 
buckling  mode  (for  example,  externally  stiffened  corrugated  cylinders  or  unstiffened 
isotropic  cylinders)  or  when  the  cylinder  has  a  pure  compressive  buckling  mode  with  a 
large  number  of  circumferential  waves  (longitudinally  stiffened  cylinders).  However,  the 
results  indicate  large  differences  between  bending  and  compression  for  cylinders  which 
have  pure  compressive  modes  consisting  of  intermediate  numbers  (for  example,  6  to  15) 
of  circumferential  waves,  such  as  the  internally  stiffened  corrugated  cylinders  or  the 
ring-  and  stringer- stiffened  isotropic  cylinders. 

Also  from  table  I  it  may  be  noted  that  the  size  of  the  determinant  needed  for  con¬ 
vergence  will  depend  upon  the  number  of  waves  in  .the  circumferential  pure  compressive 
buckling  mode;  that  is,  cylinders  buckling  with  a  small  number  of  circumferential  waves 
in  compression  will  require  fewer  terms  to  approximate  the  buckle  shape  in  bending  and 
a  smaller  determinant  order  will  be  required  for  convergence  of  the  buckling  load. 

CONCLUDING  REMARKS 

A  stability  equation  for  eccentrically  stiffened  orthotropic  cylinders  subject  to  a 
loading  of  pure  bending  or  to  combination  bending  and  compression  has  been  derived  by 
means  of  the  Galerkin  method  for  simple  support  boundary  conditions.  Sample  calcula¬ 
tions  have  been  shown  for  three  types  of  contemporary  stiffened  cylinders  to  illustrate  the 
importance  of  bending  loads,  eccentricity  effects  under  these  loads,  and  interaction 
between  bending  and  compression.  The  calculations  presented  show  that  the  difference 
in  predicted  buckling  load  between  a  cylinder  loaded  in  pure  bending  or  pure  compression 
may  be  considered  small  when  the  pure  compressive  buckling  mode  is  axisymmetric  or 
has  a  large  number  of  circumferential  waves.  The  buckling  analysis  of  cylinders,  which 
buckle  in  compressive  buckling  modes  with  a  moderate  number  of  circumferential  waves 
and  which  are  to  be  loaded  in  bending,  require  the  use  of  an  analysis  which  accounts  for 
the  bending  load.  The  calculations  also  show  substantial  eccentricity  effects  in  bending 
comparable  with  those  obtained  for  cylinders  in  compression. 

Langley  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Langley  Station,  Hampton,  Va.,  January  14,  1966. 
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TABLE  I.-  BUCKLING  CALCULATIONS  FOR  STIFFENED  CYLINDERS 
(a)  Ring- stiffened  corrugated  cylinders  (b)  Ring-  and  stringer-stiffened  cylinders  (c)  Longitudinally  stiffened  cylinders 


l/R 

Nx,/Et 

m 

Minimum 

determinant* 

Nx,/Nxc 

Stiffeners,  external 

0.15 

0.007133 

1 

12 

1.013 

.20 

.004810 

1 

12 

1.015 

.25 

.003681 

1 

26 

1.019 

.30 

.002866 

1 

32 

1.018 

.35 

.002291 

1 

34 

1.022 

.40 

.001883 

1 

36 

1.022 

.45 

.001596 

1 

34 

1.027 

.50 

.001379 

1 

34 

1.029 

Stiffeners,  internal 

0.15 

o 

o 

o 

05 

CXI 

00 

00 

rr 

30 

1.021 

.20 

.003909 

1 

28 

1.027 

.25 

.002522 

1 

28 

1.032 

.30 

.001769 

1 

28 

1.037 

.35 

.001315 

1 

28 

1.047 

.40 

.001022 

1 

26 

1.048 

.45 

.000821 

1 

26 

1.053 

.50 

.000679 

1 

L 

26 

1.058 

l/R  Nx^^Et  m 

Minimum 

determinant* 

Stringers,  external;  rings,  external 

0.05  0.006239  5 

10 

1.045 

.10  .005318  4 

12 

1.038 

.15  .004972  4 

12 

1.050 

.20  .004788  4 

12 

1.072 

.25  .004666  3 

14 

1.089 

Stringers,  external;  rings,  internal 

0.05  0.004752  4 

14 

1.156 

.10  .004335  3 

14 

1.130 

.15  .004144  3 

14 

1.114 

.20  .004025  3 

16 

1.109 

.25  .003938  3 

16 

1.101 

Stringers,  internal;  rings,  internal 

0.05  0.004470  3 

14 

1.187 

.10  .003949  3 

14 

1.151 

.15  .003669  3 

14 

1.159 

,20  .003476  3 

16 

1.163 

.25  .003319  2 

16 

1.194 

l/R 

ra 

Minimum 

determinant* 

Rings,  external 

0.15 

0.005514  2 

6 

1.089 

.20 

.004679  2 

6 

1.091 

.25 

.004171  2 

6 

1.090 

.30 

.003828  2 

6 

1.088 

.35 

.003581  2 

6 

1.086 

.40 

.003391  2 

7 

1.083 

.45 

.003245  2 

7 

1.080 

.50 

.003128  2 

7 

1.079 

Ring 

s,  internal 

0.15 

0.002458  1 

10 

1.396 

.20 

.002084  1 

10 

1.397 

.25 

.001839  1 

11 

1.389 

.30 

.001666  1 

11 

1.376 

.35 

.001536  1 

11 

1.376 

.40 

.001442  1 

11 

1.352 

.45 

.001353  1 

11 

1.341 

.50 

.001286  1 

11 

1.329 

_ 

^Minimum  order  of  determinants  required  for  0.5-percent  error  in  buckling  load. 
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R  =  200  (508) 


p  =4.00(10.16) 

9  =600 

t  =0.16(0.41) 
c 


0.10 

(0.25) 


Ring  section 


(a)  Ring-stiffened  corrugated  cylinders;  a  -  200  in.  (508  cm). 


LTL _ 3 

2.50  _ 
(6.35) 

. _ 3 

0.10  J 

1 

(0.25)  ' 

R  =  200  (508) 


Stringer 


0.90 

(2.29) 

0.10 

(0.25) 


T 

2.00 

(5.08) 

0.90 

(2.29) 


Stringer  section 


(b)  Ring-  and  stringer-stiffened  cylinders;  a  =  200  in.  (508  cm).  (See  fig.  2(a)  for  ring  section.) 

Figure  2.-  Dimensions  of  stiffened  cylinder  walls.  Dimensions  are  in  inches.  (Parenthetical  dimensions  in  cm.) 
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Maximum  pure  bending 


Figure  3.-  General  instability  predictions  for  ring-stiffened  corrugated  cylinders 


Figure  4.-  General  instability  predictions  for  ring-  and  stringer-stiffened  cylinders. 


Maximum  pure  bending 


Figure  5.-  Panel  instability  predictions  for  longitudinally  stiffened  cylinders. 


.8 


1.0 


Figure  6.-  Interaction  curve  for  ring-stiffened  corrugated  cylinder  with  l/R  =  0.2. 
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Pure  compression  (ref.  3) 


Figure  9.-  Predicted  mode  shapes  for  ring-  and  stringer-stiffened  cylinder  with  l/R  =  0.15. 
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